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Abstract 

Neutrinoless double beta decay is the most powerful tool to probe not only for Majo- 
rana neutrino masses but for lepton number violating physics in general. We discuss 
^ ■ relations between lepton number violation, double beta decay and neutrino mass, 

review a general Lorentz invariant parametrization of the double beta decay rate, 
'kj^- highlight a number of different new physics models showing how different mecha- 

: "^^"^^ ^""^^^ possibilities to discriminate 

■ and test these models and mechanisms in complementary experiments. 
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1 Introduction 

The search for neutrinoless double beta decay (Oz//3/3) - the simultaneous transformation of 
two neutrons into two protons, two electrons and nothing else - is the most sensitive tool 
for probing Majorana neutrino masses (see the contribution by Rodejohann pQ). However, 
while this so-called mass mechanism is certainly the most prominent realization of the decay, 
and while an uncontroversial detection of neutrinoless double beta decay will inevitably 
guarantee that neutrinos are Majorana particles, Majorana neutrino masses are not the 
only element of beyond Standard Model physics which can induce double beta. 
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In this review we discuss mechanisms of neutrinoless double beta decay where the 
lepton number violation (LNV), necessary for the decay, does not directly originate from 
Majorana neutrino masses but rather due to lepton number violating masses or couplings 
of new particles appearing in various possible extensions of the Standard Model. While 
the same couplings will also induce Majorana neutrino masses, due to the Schechter-Valle 
black box theorem OE], in these cases the double beta decay half life will not yield any 
direct information about the neutrino mass. 

We start our review with a general consideration of black box contributions to neu- 
trinoless double beta decay and neutrino masses. Next, we discuss a general framework 
which allows to parametrize and analyze any single contribution to neutrinoless double 
beta decay allowed by Lorentz invariance. While the neutrino mass limit is based on the 
well-known mechanism exchanging a massive Majorana neutrino between two standard 
model {V — A) vertices, the effective vertices appearing in the new contributions involve 
non-standard currents such as scalar, pseudoscalar and tensor currents. Moreover, besides 
contributions with a light neutrino being exchanged between two separated vertices, the 
so-called long-range part, additional contributions from short range mechanisms are possi- 
ble, where exchanged particles are all much heavier than the typical length scale of nuclear 
separation, such as in supersjTumetry (SUSY) without i?-parity. Then we turn to several 
concrete models such as left-right symmetry |11[5], i?-parity violating SUSY |6HT0] and 
leptoquarks p!T|fT2] . Finally we discuss the prospects to discriminate different mechanisms. 

2 Black Box Theorem 

The observation of neutrinoless double beta decay demonstrates that lepton number is 
violated. Lepton number violation implies that neutrinos have to be Majorana particles. 
That the two are inseparably connected can be proven by what is known as the black box 
theorem [2l[3l[T3]. Graphically the theorem can be depicted as shown in Fig. [1] (left): If 
double beta decay has been seen, a Majorana neutrino mass term is generated at higher loop 
order, even if the underlying particle physics model does not contain a tree-level neutrino 
mass. 

One might wonder, if it is possible to circumvent the connection between Majorana 
neutrino mass and double beta decay by tuning, for example, in a given model the tree-level 
and 1-loop contributions to the neutrino mass in such a way that the resulting observable 
neutrino mass is too small to be detected and this is of course possible. However, a tuning 
order-by-order in perturbation theory requires a stabilizing symmetry for it to be natural 
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Figure 1: Black box theorem graphically. To the left the original non-supersymmetric 
"butterfly" diagram, to the right the supersymmetric theorem. 

and the black box theorem [2] states that no such symmetry can exist, as has been shown 
formally in |13] . 

A word of caution. The black box theorem has often been mis-presented in the 
literature. It does not guarantee by any means that the mass mechanism of neutrinoless 
double beta decay is dominant. This can be seen already by looking at the "butterfly" 
diagram in Fig. [1] (left): If this diagram is the only contribution to the neutrino mass 
the resulting rrii, is nearly infinitesimal, due to the 4-loop suppression factor. A recent 
calculation of the diagram in [H] quotes a neutrino mass of the order of mj, ~ 10^^"^ eV 
from current double beta decay limits. However, this does not imply that the neutrino mass 
has always to be this tiny in order for non-mass-mechanism contributions to double beta 
decay to dominate. In fact it is easy to find examples for both kinds of particles physics 
models, those that do give mass mechanism dominance and those that do not. The classic 
example of the former is the seesaw mechanism, which leaves (in a non-supersymmetric 
world) double beta decay dominated by neutrino mass as the only experimental signature. 
In the opposite class falls supersymmetry with (trilinear) i?-parity violation, see below. 

One can also prove a supersymmetric version of the black box theorem [T5HT7]. This 
is depicted graphically in Fig. [1] (right). In a supersymmetric theory the scalar partner of 
the ordinary neutrino is a complex field. Once lepton number is violated, this complex field 
splits into its real and imaginary components with a non-zero mass difference ml^. This 
mass splitting generates lepton number violating effects such as sneutrino to anti-sneutrino 
oscillations p^[T9] . In analogy with the ordinary black box theorem one can show that such 
a LNV mass splitting in the sneutrino sector leads to double beta decay and observation 
of double beta decay implies that is different from zero. At the same time, Majorana 
neutrino masses and m^^ are also always connected and, as shown in [T51 - IT7] , the existence 
of one implies the existence of the other (if SUSY indeed is realized). 
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Finally, we mention that the original version of the black box theorem [2l|3l[l3] con- 
siders only the first generation of leptons. Oscillation experiments have shown, however, 
that flavour violation exists in the neutrino sector. It is then possible to prove an extended 
black box theorem |20j, which takes also into account flavour. In essence, this theorem 
states that under the assumption that all three light neutrinos are Majorana particles, cur- 
rent oscillation data requires the 0^(3(3 decay observable mge, see eq. below, must be 
different from zero, since no symmetry can exist which guarantees this entry in the neu- 
trino mass matrix to vanish exactly. While academically amusing, however, the theorem 
can not predict the actual value of mee, exactly as the original version can not guarantee 
the dominance of the mass mechanism. Thus, the results of [20] do not guarantee that the 
double beta decay amplitude is observably large. 



3 Lorentz-invariant Description of Neutrinoless Dou- 
ble Beta Decay 



We continue by considering the neutrinoless double beta decay rate in a general framework, 
parametrizing the new physics contributions in terms of all effective low-energy currents 
allowed by Lorentz-invariance. This parametrization has been developed in [211|22]. Such 
an ansatz allows one to separate the nuclear physics part of double beta decay from the 
underlying particle physics model, and derive limits on arbitrary lepton number violating 
theorie^]. 

Before discussing the general decay rate, let us recall that the mass mechanism of 
double beta decay, measures (or limits) the effective Majorana neutrino mass defined as: 



Here, the sum runs over all light neutrinos j with couplings to the electron and a SM 
VT-boson. It is straightforward to show that this quantity is equal to the (ee) entry of 
the Majorana neutrino mass matrix in the basis where the charged lepton mass matrix 
is diagonal. For this reason {my) is often also denoted as rriee- This contribution to 
neutrinoless double beta decay is discussed at length in the contribution by Rodejohann [1] 
to this issue. 

While the general decay rate is independent of the underlying nuclear physics model, 
to extract quantitative limits values for nuclear matrix elements are needed. Limits dis- 
cussed below are derived using matrix elements calculated in proton-neutron (pn) QRPA. 

^For another approach based on an effective operator description compare [25] . 




(1) 
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For the isotope ''^Ge this matrix elements were already available in the literature [71 I2T1I22] . 
For other isotopes, quoted in the tables below, the numerical values are taken from [21]. 

Uncertainties in nuclear matrix elements are notoriously difficult to estimate and 
all limits derived from double beta decay suffer from these uncertainties. Unfortunately, 
despite all the efforts devoted to the improvement of the matrix element calculations, the 
latest QRPA matrix elements from the Tiibingen group f25] differ from the shell model 



results revisited 



or example, in [26] in many cases by factors of ~ (2 — 3) in case of the 
Moreover, shell model matrix elements are up to now available only for 
the mass mechanism. Thus uncertainties in other matrix elements, needed in the general 
decay rate, are even harder to estimate. However, for the long-range part of the amplitude, 
discussed in section (13.11) we believe that all matrix elements suffer from uncertainties of 
the same order as those found for the mass mechanism. 

For the short-range part of the amplitude, see section (13.21) . no other general calcula- 
tion than the one presented in [22] exists. However, [25] contains matrix elements for heavy 
neutrino exchange and for the short-range /2p-violating SUSY mechanism, which we can 
compare to [22]. One noticeable difference is that in [28l|29] it was argued that the effect 
of short range correlations had been overestimated in earlier calculations. A recalculation 
of the nuclear matrix elements in [30], using the method proposed in [28], indeed led to an 
increase of (25-40) % in the numerical values of the nuclear matrix elements. The latest 
calculation [25] also has short range matrix elements which are larger than those in [22] by 
similar factors. Despite these more recent calculations we will stick to the matrix elements 
presented in [71|211|22l[2l] , since (a) no other complete calculation of matrix elements exist 
and (b) newer matrix elements in existing cases tend to be larger than those of the above 
publications, i.e. we believe that our limits are conservative. 

Currently the most stringent bounds on neutrinoless double beta decay come from 
^^Ge [31] and ^^^Xe [32]. Starting from around 2001 [33H35] a small part of the Heidelberg- 
Moscow collaboration claimed to have observed evidence for neutrinoless double beta decay, 
but this has so far not been conffimed in any other experiment. In fact, the recent publi- 
cation of the limit from ^'^^Xe [32] puts some pressure on the claim, although, due to the 
uncertainty in the nuclear matrix element calculation, ^^^Xe can not unequivocally rule it 
out yet. For definiteness we will use the limit from ^^Ge of T1/2 > 1-9 x 10^^ ys [31] and 
the recent result T1/2 > 1-6 x 10^^ ys for ^^^Xe [32] for the derivation of limits. Results for 
these two isotopes lead currently to very similar limits, see below 



^The calculation for ^^^^Xe is a notable exception. Here, the latest shell model matrix element for the 
mass mechanism agree with |27] within the error bars estimated for the QRPA calculation. 

•^For the mass mechanism, using the nuclear matrix elements from |36j . the ''^Ge limit corresponds to 
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Figure 2: Different contributions to the general double beta rate: The contributions (a) - (c) 
correspond to the long range part, the contribution (d) is the short range part (from |21j). 
(a) corresponds to the mass mechanism. 

3.1 Long— Range Part 

This subsection is essentially based on reference [21]. We consider first the long-range part 
of neutrinoless double beta decay with two vertices, which are pointlike at the Fermi scale, 
and exchange of a light neutrino in between. The general Lagrangian can be written in 
terms of effective couplings e^, which correspond to the pointlike vertices at the Fermi scale 
so that Fierz rearrangement is applicable, 

with the combinations of hadronic and leptonic Lorentz currents = uOad and jp = eOpv 
of defined helicity, respectively. The operators Oa,p are defined as 

Ov-A = - 75), Ov+A = r{i + 75), 

= (1-75), 05+P = (l + 75), (3) 

Ot, = ^[7m,7.](1 - 75), Or, = '-b„lu]{l + 75). 

The prime indicates the sum runs over all contractions allowed by Lorentz-invariance, 
except for a = (3 = {V — A). Note that all currents have been scaled relative to the 
strength of the ordinary {V — A) interaction. 

(m^) < 0.35 eV, while the ^^^Xe gives (m^) < 0.34 eV. 
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The effective Lagrangian given in eq. (Ej) represents the most general low-energy 4- 
fermion charged-current interaction allowed by Lorentz invariance. The interpretation of 
the effective couplings e^, however, depend on the specific particle physics model. Nev- 
ertheless one realizes the following general feature. Using only the SM fermion fields and 
working in the Majorana basis for the neutrinos (z/ := + u^) it is easily seen that all 
currents involving operators proportional to (1 + 75) violate lepton number by two units, 
i.e. the corresponding ej^ must also be lepton-number violating. Such LNV are easily 
found, an example is given by i?-parity violating supersymmetry treated in ref. [7t[T0|[37] 
and discussed later. 

The double beta decay amplitude is proportional to the time-ordered product of two 
effective Lagrangians (see Fig. |2]): 

T(£(i)/:(2)) = ^T{3y_A4_^3v-A4.A + ^3 pJ^Jv ^A^^A + e?J/34jV^ }• (4) 

The first term (Fig. [2] (a)) corresponds to the contribution from the Majorana neutrino 
mass, and the 3rd term (Fig. [2] (c)), which is quadratic in e can be neglected. Only the 
2nd term (Fig. [2] (b)) is phenomenologically interesting. For this term one has to consider 
two general cases: 

1) The leptonic SM {V — A) current meets a left-handed non SM current jp with /3 = 
{S — P),Tl. For this contribution the neutrino propagator is 

with the usual left- and right-handed projectors Pl/r = This expression is propor- 

tional to the unknown Majorana neutrino mass rriu < 0.5 eV, for which no lower bound 
exists. Therefore no limits on the corresponding parameters can be derived. 

2) The leptonic SM {V — A) current meets a right-handed non SM current jp with /3 = 
{S + P),(y + A),Tji. For this contribution the neutrino propagator is 

q^%^ + m^ _ q^l^, 
qf.q" - ml q^qt" - mf, 

which is proportional to the neutrino momentum. Since typically q/j^ c:^ pp — 100 MeV with 
the nuclear Fermi momentum pp, this part of the amplitude will produce stringent limits 
on corresponding e^. 

Considering only one at a time (evaluation "on axis") one can now derive con- 
straints on the effective coupling parameters from a double beta decay half life measurement 
or bound, 

[r°;f = |6^^Go,|ME|^ (7) 
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where Gok denotes the phase space factors given in [4J and \ME\ the nuclear matrix elements 
discussed below. Note that evaluating "on axis", compared to the arbitrary evaluation, 
neglects interference terms of the different contributions. Numerical values of the matrix 
elements J^j.', A^^^', A^^^", A^,^^", J^p' below, calculated in the Quasi Particle Random 
Phase Approximation (pn-QRPA), are given in Table O The matrix elements for the 
different contributions are defined as follows. 



SM meets jv+aJv+a Jv+aJv-a These combinations of currents and the corre- 
sponding matrix elements have been discussed in the literature before [ll[5l|36]. Matrix 
elements have been calculated in different papers, and we will follow [5]. 



SM meets js+pJs+p ^ind js+pJs-p Using s-wave approximation for the outgoing elec- 
trons and some assumptions according to 0,[38] one has 



_FfiO)f 1 

with the phase space factor Gqi and the Gamov- Teller and tensor matrix elements M-gt' 
and J^rp', respectively. In addition, R denotes the nuclear radius, rrie the electron mass, 
Ga ^ 1.26 and Fj^\o) = 4.41 §9\. 



SM meets JtrJtj^ and Jt^Jt^ The hadronic Tr contribution is given by 

MEg = -ai'^Mar' + aiM^', (9) 

with an effective nuclear form factor ai. For the hadronic contribution in leading order 
of the inverse proton mass {l/rrip) one finds 

ME^^ = a2Mp' - a-s (^M^" + ^M^t") , (10) 

with nuclear matrix elements Mp', Mqp" and Aip" and effective nuclear form factors a2 
and 0:3. The parameters are defined as: 

4:Ti'\0)Gvil-2mp{Gw/Gv)) 
= G^, ' 

4(2ff (O)-Tf (0))Gv. 
= G^, ' ^^'^ 

ATi'\0){Gp/GA) 
= GAR^m, ■ ^'^^ 
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Isotope 


G 


01 




G 


06 




G 


09 






6.40- 


10- 


-15 


1.43 ■ 


10" 


-12 


3.30 • 


10- 


-10 




2.82 ■ 


10" 


- 14 


4.77- 


10" 


-12 


1.32 ■ 


10" 


-9 


looMo 


4.58- 


10- 


-14 


7.09- 


10- 


-12 


1.88- 


10- 


-9 


i28Te 


1.83- 


10" 


-15 


4.94- 


10- 


-13 


7.52- 


10- 


-11 




4.44- 


10- 


-14 


6.89- 


10- 


-12 


1.55- 


10- 


-9 




4.73- 


10- 


-14 


7.28- 


10- 


-12 


1.60- 


10- 


-9 


150Nd 


2.10 ■ 


10- 


-13 


2.48 ■ 


10- 


-11 


6.45- 


10- 


-9 



Table 1: Phase space factors for the general decay rate, numerical values taken from the 
calculation of [3]. 

Here, t[^^ = 1.38, ^ = -4.54 [39j, Gp/Ga = Im^jml and {GwlGy) = - st; 

obtained from the CVC hypothesis. Numerical values for the matrix elements and limits 
on the eg are discussed in Section 13.31 



3.2 Short-Range Part 

In the short range part the effective interaction can be considered as point-like, thus the 
decay rate results from the following general Lorentz invariant LagrangiaiJ^: 

>C = ^m-\t^JJ3^t2J^-J^,3^t;J^J^3^t,J^J^,3- ^t^J^Jj^}, (14) 

with the hadronic currents of defined chirahty J = u(l ± 75)^, = ?l7^(l ± 75)^, J'^'^ = 
u|[7^,7'^](l ±75)(i and the leptonic currents j = e(l ±75)6*-^, = 67^(1 ±75)6*". In some 
of the cases the decay rate for the effective coupling depends also on the chirality of the 
currents involved. In these cases we define = e^^^, where xyz = L/R, L/R, L/R defines 
the chirality of the hadronic and leptonic currents in the order of appearance in eq. ( lTi|) . 
In the cases where it is not necessary to distinguish the different chiralities we suppress 
this additional index. 

In renormalizable theories no fundamental tensors exist. Thus the tensor currents 
have to result either from Fierz rearrangements or by integrating out heavy particles, 
when deriving the effective Lagrangian from the fundamental theory and decomposing the 
expressions obtained in terms of the Lorentz invariant bilinears used above, e.g. from 
'U7^7,^(l + 75)(i = g^uJ — iJfiu- Applying the standard nuclear theory methods based on the 
non-relativistic impulse approximation one derives the general Oz//3/3-decay half-life formula 



*Here we follow the essentially the calculations presented in [2^, see however, Table [2j 
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Isotope 




Mf' 


■Mgt" 


Mt' 


Mt" 


■M.GT,N 


Mf,n 




2.95 


-0.663 


8.78 


0.224 


1.33 


0.113 


-0.0407 




2.71 


-0.603 


7.96 


0.208 


1.26 


0.102 


-0.0360 


looMo 


3.69 


-0.876 


13.4 


0.328 


2.44 


0.129 


-0.0489 




2.26 


-0.509 


7.13 


0.193 


1.47 


0.075 


-0.0271 


i28xe 


3.70 


-0.814 


13.4 


0.331 


2.37 


0.119 


-0.0419 




3.27 


-0.720 


12.0 


0.304 


2.20 


0.105 


-0.0369 




1.83 


-0.403 


6.90 


0.165 


1.18 


0.058 


-0.0203 


150Nd 


5.39 


-1.21 


21.8 


0.642 


5.07 


0.165 


-0.0591 



Table 2: Nuclear matrix elements for Oz//3/3 decay calculated in the pn-QRPA approach. 
Results of different publications [71[2T],[22l|2l] are summarized here. Note that the isotope 
^^°Nd has a sizeable deformation, but the calculation is performed in the spherical limit. 
The numbers for ^^°Nd might therefore be an overestimation of the true values. 



in s-wave approximation 

3 



/2 J 



i=l 



Go 



5 

E 

i=4 



j=i 



i=4: 



■ (15) 



Here the phase space factors are 
Gi = Goi, G2 



09; 



G, 



rripR 



G, 



06) 



(16) 



with Gofc calculated in f4| and shown in Table [T] for completeness. The nuclear matrix 
elements in eq. f[T5|) are defined ajf] 



Ml 
Mi 



-af^MF,N, 



M, 



al'^McT, 



m 



^ {Mgt,nToiI^Mf,n} 

e 

SR 



±a^''MGT,N 



M, 



N, 



N- 



(17) 



In the last three cases contractions of hadronic currents with different chiralities lead to 
different results. The negative sign in M^ corresponds to 63^^ and ef^"" {.JvtaJvta)i the 
positive sign to e^"^^ and ef^^ {JvtaJv±a)- The sign of Mi is positive for the combinations 

{JvtaJtl/trJv-a and Jv±AJTL/TRjv+A)- For the combinations 



^LLL ^FtRL ^RLR ^LRFt 



r-LLR ^RRR ^RLL ^LRL 



5 



5 



, ej {JvtaJtl/trJv+a and Jv±AJTL/TRjv-A) it is negative. The sign of 



the matrix element A^s is negative for the left-handed leptonic current (eg^^) and positive 



^Note that the at in the long-range part and the af^ defined here are different coefficients. 
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Isotope 


\^V+A\ 

\^v-a\ 


\^V+A\ 

\^v+a\ 


\^S+P\ 


\^S+P\ 


1 ,rP| 

I^TL 1 


\cTR\ 




3.5-10-9 


6.2 - 10-^ 


1.1 - 10-* 


1.1 - 10-* 


6.7-10-1° 


1.1 - 10-9 


i36Xe 


2.8 - 10-9 


5.6-10-^ 


6.8- 10-9 


6.8-10-9 


4.8 - 10-1° 


8.1 - 10-1° 



Table 3: Limits on effective long-range B — L violating couplings. These limits are derived 
assuming only one e is different from zero at a time. With the recent limit on the half-live 
for i^^Xe [32], i^^Xe now gives limits competitive with or better than ''^Ge. 



for the right-handed one (eg^^). The numerical values of the standard nuclear matrix 
elements M.f,n and J^gt,n in eq. (fTTll . calculated in the Quasi Particle Random Phase 
Approximation (pn-QRPA), are given in Table [21 The pre-factors af^ in eq. ( JT7j) are 
defined as follows, 



V Ga ' mpTrie ' V Ga ^ mpm, 



P I I N V I ' II Lip I I VQ 

^SR_(9v_V ^SR _ T^A ^SR _ 9vFP ml 



3 \GaJ ' ' GAmpm^ ' G^ 

The finite nucleon size is taken into account in a common way [lOll^T] by introducing the 
nucleon form factors in a dipole form 

with tua = 0.85 GeV, gv = 1.0, = 1.26. The other form factor normalizations have 
been calculated in ref. [39] within the MIT bag model, = 0.48. 



3.3 General Oi/f^P Constraints 

We list all nuclear matrix elements necessary for deriving limits for the eight most important 
nuclear isotopes in Table |5J The long range NMEs for ^^Ge had been published previously 
in [21]. The long-range matrix elements for other isotopes are from [21]. Note that the 
numerical values for i^°Nd might overestimate the true NMEs, since the calculation was 
done in the spherical limit. With these matrix elements we find the limits on the different 
contributions presented in Tables [3] and ID 

A few comments on tables ([3]) and (jl]) might be in order. The Jv+aJv+a 
Jv+aJv-a P^^t of the amplitude has been considered within left-right symmetric mod- 
els [ll[5l|36]. In our notation, the limits are given as e^^^ and ey^^. In the notation of [1] 
these correspond to (A) = ej^^^ and (77) = ey'^^. Note also, that we have updated the 
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ei 




£2 


1 LLziRRz), 


1 LRz(RLz), 
F3 1 




£4 




£5! 


™Ge 


3.2 


•10-7 


1.8 • 10-9 


2.2 • 10-8 


1.4 • 10-8 


1.5 


•10-8 


1.5 


•10-7 


i36Xe 


2.6 


•10-7 


1.4 • 10-9 


1.1 • 10-8 


1.7 • 10-8 


1.2 


• 10-8 


1.2 


•10-7 



Table 4: Limits on effective short-range B — L violating couplings. These limits are derived 
assuming only one e is different from zero at a time. For 63 the contractions of hadronic 
currents with different chiralities lead to different results. With the recent limit on the 
half-live for ^^^Xe [32], ^^^Xe now gives limits competitive with or better than "^^Ge. 

limits with the half-live limits from [31] for 76Ge and for ^^^Xe from [32] ■ AH limits are 
derived "on axis" , i.e. assuming only one non-zero contribution at a time. 

4 Models of Lepton Number Violation 

In the following we discuss several prominent new physics models that incorporate lepton 
number violation and which lead to potentially observable rates for neutrinoless double 
beta decay. The list of models presented here is not intended to be exhaustive, as there 
is a large number of alternative schemes which produce interesting Oi//3/3 phenomenology 
such as scalar bilinears [42] or a scalar octet seesaw mechanism [ISjHl]. A large range of 
models is discussed in the review [12] and in the references therein. 

4.1 Left-Right Symmetry 

As a first example of a model incorporating a rich phenomenology of lepton number viola- 
tion, we will discuss the minimal Left-Right symmetric model (LRSM) which extends the 
Standard Model gauge symmetry to the group SU(2)l ® SU(2)ij ® U(1)b-l [461 - 149] . Right- 
handed neutrinos are a necessary ingredient to realize this extended symmetry and are part 
of an SU(2)i^ doublet. In the LRSM, a generation of leptons is assigned to the multiplets 
Li = iyi^li) with the quantum numbers = (1/2,0, —1) and Ql^ = (0, 1/2, —1) under 
SU(2)i ® SU(2)/j ® \]{\)b-l- The Higgs sector contains a bidoublet and two triplets A/, 
and A/j. The VEV vr of A/j breaks SU(2)/j ® U(1)b-l to U(l)y and generates masses for 
the right-handed Wr and Zr gauge bosons, and the heavy neutrinos. Since right-handed 
currents and particles have not been observed, vr has to be sufficiently large. The neutral 
part of the bidoublet acquires a VEV v at the electroweak scale thereby breaking the SM 
symmetry. The LRSM can accommodate a general 6x6 neutrino mass matrix in the basis 
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Figure 3: Oz//3/3 in the LRSM: Light (left) and heavy (right) neutrino exchange. 



M={^;^ ,"1, (20) 



Ml Md 
Ml Mr 

with Majorana and Dirac mass entries of the order M^ ^ VmVl, Mr ^ VmVr and Md = 
Udv. Here yM,D are Yukawa couphngs and vl is the VEV of the left Higgs triplet, which 
together with the other vacuum expectation values satisfies vlVr = v"^. The mass matrix 
( l20l) is diagonalized by a mixing matrix of the form 

with the 3x3 block matrices U and V describing the mixing among the light and heavy 
neutrinos, respectively, whereas W yields left-right mixing between the light and heavy 
states. 



4.1.1 Neutrinoless Double Beta Decay 

In the LRSM, several mechanisms can contribute to Oz//3/3 as shown in Figs. [3l [Hand [51 The 
contributions in Figs. [3] and H] are of the same diagramatical form with the exchange of either 
light or heavy neutrinos as well as light and heavy W bosons. Diagram [3] (left) describes the 
standard mechanism of light neutrino exchange, with the effective mass mee = \ Xli U^i'^i^i I > 
saturating current experimental bounds if the light neutrinos are degenerate at a mass scale 
m^-^ ^ rriee ^ 0.3 — 0.6 eV. Correspondingly, diagram [3] (right) describes the exchange of 
heavy right-handed neutrinos. In the classification of Section [HI this is a realization of 
the short-range operator with the effective coupling ef^^. Assuming manifest left-right 
symmetry, i.e. qr = qr, in terms of the LRSM model parameters it is given by 

er = EK^^^, (22) 
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Figure 4: Oz//3/3 in the LRSM: Light neutrino exchange not proportional to the neutrino 
mass, with a left-handed (left) and right-handed (right) hadronic current. 




Figure 5: Ou/S/S in the LRSM: Doubly-charged Higgs Triplet exchange. 

and searches for Ou/3(3 yield the limit jef^"^! < 1.1 ■ 10"® (cf. Table |4]). 

Because of the presence of right-handed currents in the LRSM, light neutrino exchange 
does not necessarily require a chirality violating mass insertion. As a consequence, the 
outgoing electrons can have opposite chirality, and the contributions are suppressed either 
by the heaviness of or the smallness of the mixing angle ( of the W bosons as shown 
in Fig. HJ The coupling parameters of the corresponding effective long-range operators can 
be written as 

3 2 3 

4Xi = Yl UeiW.i'^, elti = J2 tan C, (23) 

with the current experimental limits |e|/+^| < 5.6- 10"'^ and |e|/l^| < 2.8- 10~^, respectively 
(cf. Table [3]). Both cases are necessarily suppressed by the left-right neutrino mixing 
Md/Mn ~ yjmyjmf^ (the latter expression is valid for a dominant type-I seesaw mass 
mechanism [50]) between light and heavy neutrinos. 

Finally, Fig. [5] describes the exchange of a right-handed doubly-charged triplet Higgs 
A/{, which has the same effective operator structure as heavy neutrino exchange. The 
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Figure 6: Production and decay of a heavy right-handed neutrino with dilepton signature at 
hadron colhders (left). Comparison of LNV event rates at the LHC and in Oz//3/3 experiments 
(right, from [51]). The sohd blue contours give the number of e^e^ + 2j events at the LHC 
with 14 TeV and £ = 30 fb~^. The dashed blue contours correspond to signal significances 
of 5cr and 90%. The shaded green area denotes the parameter space excluded by 0^(3(3 
at T^^^l^ ^ 2 X 10^^ years, assuming dominant doubly- charged Higgs or heavy neutrino 
exchange. The green dashed contours show the sensitivity of future Oz//3/3 experiments at 
j^Qupp ~ 2 X 10^^ years. The red shaded area is excluded by LHC searches [52] . 

effective short-range coupling strength is here given as 



currently limited to < 1.1 x 10^*^ (cf. Table |4]). The heavy neutrino masses m^r. 

appear since in Figure [5]the coupling of the Higgs triplet to the gauge boson is proportional 
to vr and the electron vertex is of Yukawa strength {yM)ee- 

4.1.2 Lepton Number Violation at the LHC 

In [51] , the potential to discover number violating dilepton signals pp — )■ Wr — )■ e^/x^'^ + 
2 jets via a heavy right-handed neutrino at the LHC was determined, cf. Fig. [6] (left). 
Fig. ini (right) compares the LHC event rates for such lepton number violation processes 
with the sensitivity of Qvfifi experiments. The green regions and green dashed contours 
represent the excluded areas from Qv^fi searches using nominal values for the current limit 




(24) 
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(^rpovpii > 2 X 10^^ years) and the future sensitivity (T°^^^ ^ 2 x 10^® years). In this analysis 
it was assumed that left-right mixing is negligible and Oz//3/3 is either dominated by heavy 
neutrino or Higgs triplet exchange. As the contribution from the standard light neutrino 
exchange is always present, the these results correspond to a scenario with a small effective 
mass niee- Fig. [6] (right) gives an example of the possible interplay between searches for 
lepton number violation at a high energy collider and in Qufifi experiments. 

4.2 i?— Parity Violating Supersymmetry 

The MSSM (minimal supersymmetric extension of the standard model) assumes a discrete 
Z2 symmetry, called i?-=parity (-Rp), exists. This symmetry guarantees the lightest super- 
symmetric particle to be stable, thus the MSSM with Rp offers a dark matter candidate. 
However, from a theoretical perspective the MSSM does not offer any explanation as to why 
Rp is conserved. Rather it is an ad hoc symmetry to avoid a phenomenological disaster. 
Consider the Rp violating terms 

W$.p = \jkLiLjEk + X[j^LiQjDk + eiLiHu + X'ljjJJiDjDk, (25) 

where indices 2, j, k label generations. The first three terms violate lepton number, while the 
last one violates baryon number. In the presence of both types of terms the proton decays 
at a rate which is many orders of magnitude above the experimental bound. However, any 
discrete symmetry which eliminates either the baryon or the lepton number violating terms 
is phenomenologically acceptable [53]. In fact, since the lepton number violating terms in 
eq. fl25l) generate Majorana neutrino masses, a small amount of Rp violation could actually 
explain the observed neutrino oscillation data [51] . 

4.2.1 Neutrinoless Double Beta Decay 

In the case that Rp is broken, Oi//9/3 decay can occur through Feynman graphs involving 
the exchange of superpartners as well as -^p-couplings A' [SI [71 [HI [TD]. The short-range 
contribution has been discussed in 16, 7|. Here, we take the opportunity to correct some 
errors in the original publication [7] , which leads to a slight change in the numerical values of 
the published matrix elements: (i) In the coupling constants «^/y of [7], eqs. (57)- (60), one 
should replace Tlf"* — )■ Tlf' . (ii) In Section V of [7j five different nuclear matrix elements are 
defined. The tensor matrix element should be written as A^^' ~ (o"i ■ fij)(o"j " fij) — l/SfXi ■ ay 
And, finally (iii) there was a numerical error in the code, which converted Table II to 
Table III in [7j. The corrected values for the nuclear matrix element M.q are given in 
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^^^Ge 82se loo^o ne^d iss^e isoje is^Xe ^^^Nd 



A) Mq 222 199 265 

B) Mg 281 251 331 



190 



150 



242 



301 



266 



214 



145 



116 



432 



353 



Table 5: Nuclear matrix elements for short-range SUSY Ou/S/S decay. Shown are Aiq for 



errors discussed in the text. 

Table [51 As in the original paper, —Aiq is given, to account for a relative sign in the 
neutrino mass mechanism with respect to the definitions used in the formalism of Doi, 
Kotani and Takasugi Note that the Rp violating SUSY mechanism depends on a 
combination of different short-range matrix elements discussed in section (13. 2p , since in the 
amplitude both JJj and J^^ J^vj currents appear. 

Ref. [7] used the half-live limit from the Heidelberg-Moscow collaboration available at 
that time [55] , which was later superseded by the more stringent value of |31] . The change 
in the matrix element and the update in the half-life combined leads to a slightly more 
stringent limit on A'^^^^ given by 



for mj^ = rrtu^. For comparison, from ^^^Xe data [32] one gets currently A'^^ii < 2.7 ■ 10~^. 
Note, that recently [25] published sets of matrix elements for Rp violating double beta 
decay for the pion-exchange mechanism, which are, depending on the choice of nucleon- 
nucleon interaction, model space and value of Ga-, between a factor of x ~ (2 — 3) larger 
than the matrix elements of Table O These would lead to limits on A'^^ii which are stronger 
by a corresponding factor of y/x. 

0u/3(3 decay is not only sensitive to A'^^^^^. Taking into account the fact that the SUSY 
partners of the left- and right-handed quark states can mix with each other, new diagrams 
appear in which the neutrino-mediated double beta decay is accompanied by SUSY ex- 
change in the vertices [814T0]. A calculation of previously neglected tensor contributions 
to the decay rate allows to derive improved limits on different combinations of A' [TO] . 
Assuming the supersymmetric mass parameters of order 100 GeV, the half life limit of the 
Heidelberg-Moscow Experiment implies: A'^^^gA'-^g^ < 3 ■ 10~^, A'^]^2^'i2i — ^ ' 10~^. 



the two sets of input values of coefficients for the a^*^ of Table 1 in [7], corrected for the 




(26) 
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4.2.2 Lepton Number Violation at the LHC 

Similar to the situation in left-right symmetric models, also in i?-parity violating SUSY 
scenarios interesting and complementary information can be obtained from the neutrinoless 
double beta decay analogue at the LHC, namely resonant single selectron production with 
two like sign electrons in the final state [56|[57]. The color and spin-averaged parton total 
cross section of a single slepton production is given by |58] 

where s is the partonic center of mass energy, mj is the mass of the resonant slepton, 
and finite width effects have been neglected. Considering effects from parton distribution 
functions, to a good approximation the total cross section scales like a{pp — ?■ /) oc l^'ml'^ /itij 
with the slepton mass in the parameter region of interest. 

Thus for small slepton masses, the stringent bound from Oz//3/3 decay makes this 
process unobservable at the LHC. However, the bound on A'^^^^^ originating from the non- 
observation of Oz//3/3 decay scales with the slepton mass like a < cAgjjgy where c is a 
constant, so that for higher values of the SUSY masses, larger cross-sections may be allowed 
as much a larger A'^^^^^^ is no longer excluded. It is this possibility that can be exploited at 
the LHC. 

In [SBIIFT] it has been shown that much of the parameter space allowed by Oi//3/3 decay 
in simple models of supersymmetry breaking actually predicts observable single slepton 
production at the LHC (compare Fig. [7]). Moreover, if the next generation of experiments 
observe Oz//3/3 decay, the LHC has a very good chance of observing single slepton production 
with only 10 fb~^ of integrated luminosity, assuming that Oz//3/3 decay is induced by a A'^^^^^ 
coupling. On the other hand, non-observation of single slepton production could then 
discriminate against the X'm mechanism. In general, both Majorana neutrino masses and 
X'lii could contribute simultaneously and non-negligibly to 0^(3/3 decay. In this case detailed 
LHC measurements of the kinematics in single slepton production could constrain the SUSY 
parameters, and the total cross-section could then give information about the size of |A']^^]^|. 
In principle and depending on nuclear matrix element uncertainties the LHC information 
could be combined to predict an associated inverse double beta decay half life coming from 
X'lii, which could be compared with the experimental half life measurement in order to see 
if additional contributions were necessary. 
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Figure 7: Production and decay of a single selectron in i?-parity violating models (left, 
from [56]) . Region in the mSUGRA parameter space (vertical axis: M1/2 in GeV, horizontal 
axis: Mq in GeV) in which single slepton production may be observed at the LHC for 
tan/3 = 10, Aq = and 10fb~^ of integrated luminosity. In the top left-hand black triangle, 
the stau is the LSP, a case not discussed in [56]. The black region at the bottom is ruled out 
by direct search constraints. The labeled contours are taken from Ref. [59], and indicate 
the search reach given by the corresponding value of X'ui- The white, dark-shaded and 
light-shaded regions demonstrate that observation of single slepton production at the 5cr 
level would imply T°/f ^ < 1.9- lO^^yrs, 100 > T°/f ^/lO^^yrs > 1.9 and r°f ^ > 1 x lO^^rs, 
respectively (right, from [56j). 

4.3 Leptoquarks 

Leptoquarks (LQs) are hypothetical scalar or vector particles coupling to both leptons 
and quarks. They appear most prominently in grand unified theories, but also in extended 
Technicolor or Compositeness models. LQs which conserve baryon number can be relatively 
light [60j, possibly within reach of accelerator experiments. Also low-energy precision 
measurements can give limits on LQ properties, for a detailed list on constraints from non- 
accelerator searches see, for example [HI] and [62] • The mixing of different LQ multiplets by 
a possible leptoquark-Higgs coupling [11] can lead to a contribution to Oz//3/3 decay, if these 
couplings violate lepton number [12]. Diagrams involving LQs and standard model weak 
current interactions can be generated, see Fig. [8l These diagrams are of the long range 
type and due to the chirality violating LQ interaction gain a ^ -enhancement in the double 
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Figure 8: Leptoquark diagrams for neutrinoless double beta decay. 

beta decay amplitude. Combined with a lower limit on the Oi//3/3 decay half-life bounds on 
effective couplings can be derived [12]. Assuming only one lepton number violating AL = 2 
LQ-Higgs coupling unequal to zero and the leptoquark masses not too different, one can 
derive from this limit a bound on the LQ-Higgs coupling, 

iLQ-Higgs = few ■ 10-^ (28) 

for LQ masses of the order of C(200 GeV). Such lepton number violating LQ-Higgs cou- 
plings also lead to non-zero neutrino masses at the 1-loop level [63] . 

4.4 Extra Dimensions 

Theories with large compact extra dimensions of TeV size [6H472] have enriched dramat- 
ically the perspectives of the search for physics beyond the Standard Model. Among the 
possible higher-dimensional realizations, sterile neutrinos propagating in such extra dimen- 
sions [731476] may provide interesting alternatives for generating the observed light neutrino 
masses. Conversely, detailed experimental studies of neutrino properties may even shed 
light on the geometry and shape of the new dimensions. 

The minimal higher- dimensional framework of lepton number violation considers a 
5-dimensional theory compactified on a /Z2 orbifold, in which only one 5-dimensional 
(bulk) sterile neutrino is added to the field content of the SM [77]. In this minimal model, 
the SM fields are localized on a 3-|-l-dimensional Minkowski subspace, also termed 3-brane. 
This model naturally generates small neutrino masses via a higher-dimensional version of 
the seesaw mechanism [73]. With respect to neutrinoless double beta decay an interesting 
feature of such extra-dimensional models is that the excitations of the sterile neutrino in the 
compact extra dimensions, a so-called Kaluza-Klein tower of states, contributes to the decay 
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Figure 9: The SM matter is localized on a 3-brane, while the sterile singlet neutrinos 
are allowed to propagate in the bulk. This framework naturally generates small neutrino 
masses. 

rate. As the masses of the exchanged Kaluza-Klein excitations range from small masses 
giving rise to long range contributions over the 100 MeV region up to large masses with 
short range contributions, such scenarios constitute a special case which is not described 
by the effective operator parametrization in Section 3. 

Other studies on neutrinoless double beta decay were performed within the context of 
higher- dimensional models that assume a shining mechanism from a distant brane [78] and 
of theories with wrapped geometric space [79]. In the first model Oz//3/3 decay is accompanied 
with the emission of Majorons [TS], while the model presented in [TU] predicts double beta 
decay at a rate being too small to be observable signal in running experiments. Another 
approach motivated by [71] which does not consider the effect of Kaluza-Klein states is 
described in [80]. 

The minimal higher-dimensional scenario in [77] assumes that singlet neutrinos being 
neutral under the SU(2)i(8)U(l)y gauge group can freely propagate in a higher-dimensional 
space of [1 + (3 + 6)] dimensions, the so-called bulk, whereas all SM particles are localized 
on the (3 + l)-dimensional brane. If the brane were located at the one of the two orbifold 
fixed points, the lepton number violating operators would be absent as a consequence of 
the Z2 discrete symmetry. However, if the brane is shifted by an amount a 7^ 0, these 
operators are no longer absent and this breaking of lepton number can lead to observable 
effects in neutrinoless double beta decay experiments. 

One major problem of such higher-dimensional theories is their generic prediction of 
a KK neutrino spectrum of approximately degenerate states with opposite CP parities that 
lead to extremely suppressed values for the effective Majorana- neutrino mass (mj^). How- 
ever, the KK neutrinos can couple to the W bosons with unequal strength, thus avoiding 
the CP-parity cancellations in the Oz//3/3-decay amplitude. The brane- shifting parameter 
a can then be determined from the requirement that the effective Majorana mass (m,^) 
is in the observable range. To achieve this, 1/a has to be constrained to be larger than 
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the typical Fermi nuclear momentum qp = 100 MeV and much smaller than the quantum 
gravity scale Mp, or equivalently 1/Mp <^ a < I/^f- 

The masses of the Kaluza-Klein states are obtained then by diagonalizing the infinitely 
dimensional Kaluza-Klein mass matrix and result approximately as 

TL 

m(n) ~ ;^ + ^ ' (29) 

where n is the index denoting the Kaluza-Klein excitation, R is the radius of the extra 
dimension and e is the smallest diagonal entry in the neutrino mass matrix. The mixing- 
matrix elements Bev and n follow then as [73] : 



1 , 2 2D2 , ^^'^ - (n^p\2 ' ^'^^ ' 

1 + TT^m^R^ ^ 1 1? + ^ J 

where a is the brane shift parameter and (ph is a function of a, R and the Yukawa couplings. 
The Oz//3/3-decay amplitude Tou/b/b results as 

(tti ) 

= ^^McTFim,), (31) 

where the combination of nuclear matrix elements AIgtf = A^gt — A^f sensitively de- 
pends on the mass of the exchanged KK neutrino. As soon as the exchanged KK- neutrino 
mass is comparable or larger than the characteristic Fermi nuclear momentum qp ^ 
100 MeV, the nuclear matrix element A^gtf decreases as The general expression 

for the effective Majorana-neutrino mass (m,^) in eq. ( 13T|) is given by 

^ oo 

^^•'^ = T:i 7 ^ Yl Kn^in) McTFim^n)) " McTpim^) ■ (32) 

MGTF[m,} L J 

Here the first term describes the genuine higher-dimensional effect of KK-neutrino ex- 
changes, while the second term is the standard contribution of the light neutrino u. The 
dependence of the nuclear matrix element A^gtf on the KK-neutrino masses m(„) here 
has been included in the double beta observable (m^). This leads to predictions for (m,^) 
that depend on the double beta emitter isotope used in the experiment. However, the 
difference in the predictions is typically too small to work as a smoking gun for the extra- 
dimensional mechanism of Oz//3/3 decay. A particularly interesting property of this model 
is that the effective Majorana-neutrino mass (m^) is not bounded from above by the mass 
eigenvalues of the light neutrinos: It can be close to the experimental limit even for an 
infinitesimal lightest neutrino mass m,y^ which constitutes a rather unique feature of such 
higher-dimensional brane-shifted scenarios. 
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4.5 Majorons 



In many theories beyond the Standard model in addition to the neutrinoless double beta 
decay with two electrons in the final state and nothing else, lepton number violating decay 
modes can arise where new scalar [811483] or even vector particles [H] are emitted as well: 

2n^2p + 2e- + (j),2n^2p + 2e" + 20. (33) 

Majorons have been originally introduced as Nambu-Goldstone bosons being responsible 
for breaking a global lepton symmetry and generating neutrino Majorana masses [851186] . 
However, since such classical Majorons require severe fine-tuning in order to respect the 
bounds on neutrino masses and at the same time induce an observable rate for neutrinoless 
double beta decay, several new models have been proposed in which the term Majoron 
refers more broadly to a light or massless boson with couplings to neutrinos. In comparison 
to neutrinoless double beta decay with two electrons in the final state only, these decay 
modes lead to continuous spectra for the emitted electrons which can be discriminated 
from the Standard Model 2i//3/3 decay by fitting the spectral shape. While expected decay 
rates are rather small due to suppressed nuclear matrix elements [87] experimental bounds 
from the Heidelberg-Moscow experiment using ''^Ge on various models have been derived 
in [HH]. More recently, the NEMO-3 [89] and KAMLAND-Zen [90] collaborations have 
published improved bounds on these modes, constraining the Majoron- neutrino coupling 
constant from measurements in ^^''Mo, ®^Se and ^^^Xe e.g. for classical models to (gee) < 
(0.4 - 1.9) ■ 10-^ (gee) < (0.66 - 1.7) ■ 10-^ and (gee) < (0.8 - 1.6) • 10-^ respectively, 
depending on the value of the nuclear matrix element used. For classical Majorons it has 
been shown in [91] that by combining bounds from next-generation double beta experiments 
with constraints derived from supernova energy release arguments [92] neutrino-Majoron 
couplings could be constrained down to a level of 10~^. 



4.6 Non-Standard Neutrinos 

Non-standard neutrinos can contribute to neutrinoless double beta decay either as new 
flavours in addition to three partners of the charged leptons included in the Standard 
Model, or via non-standard properties or interactions. 

Additional flavours have to be either sterile or heavier than half the mass of th Z boson 
in order to comply with the width of the Z boson. They will contribute to neutrinoless 
double beta decay via mixing with the electron flavour neutrino (compare the discussion 
in the contribution of Rodejohann). An update on a heavy 4th generation's neutrino to 
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neutrinoless double beta decay is given in [93] . It implies that such heavy active neutrinos 
have to be pseudo-Dirac neutrinos with only a tiny amount of lepton number violation. 
The amount of lepton number violation is further constrained due to wash-out effects such 
weak scale pseudo-Dirac neutrinos will have on any pre-existing baryon asymmetry |94j . 

In principle the neutrino can also be a composite object. In this case a bound on the 
compositeness scale can be obtained from the neutrinoless double beta decay half life limit 
evaluated in the mass mechanism [951 - I97] . 

In |98j it has been discussed that observable double beta decay rates in extra-di- 
mensional spacetimes could be triggered by the lepton number violation induced by virtual 
black holes violating the associated global symmetry. A somewhat related mechanism arises 
in theories with a saturated black hole bound on a large number of species. Such theories 
have been advocated recently as a possible solution to the hierarchy problem and as an 
explanation of the smallness of neutrino masses [99]. Then the violation of lepton number 
can create a potential phenomenological problem of such N-copy extensions of the Standard 
Model as again due to the lowered quantum gravity scale black holes may induce TeV scale 
LNV operators generating unacceptably large rates of e.g. neutrinoless double beta decay. 
It has been shown, however, that this does not happen in such a scenario due to a specific 
compensation mechanism between contributions of different Majorana neutrino states to 
these processes. As a result rates of LNV processes are extremely small and far beyond 
experimental reach, at least for the left-handed neutrino states |100j . 

Lorentz invariance and the equivalence principle are the most important pillars of 
special and general relativity. However, certain versions of string theories allow for or 
even predict the violation of these laws. Often a violation of Lorentz invariance (VLI) 
implies that different particles can have characteristic maximal attainable velocities. The 
difference of the velocities 6v then parametrizes the size of VLI. Similarly, the corresponding 
observable describing violations of the equivalence principle (VEP) is the difference of 
characteristic couplings 5g to the gravitational potential (p. While previous studies of 
neutrino oscillations provide very restrictive bounds in the region of large mixing |101j 
Oi//3/3 decay yields in certain models of VLI and VEP a bound also in the region of zero 
mixing being not accessible to neutrino oscillation experiments |102j : 5v < 3.3 ■ 10~^^, 
(f) ■ 6g < 3.3 • 10^^^. Lorentz invariance is also closely related to CPT invariance. In 
CPT violating models the conventional notion of Majorana neutrinos being their own anti- 
particles does not apply anymore. Majorana masses and double beta decay in such scenarios 
are discussed in |103] . 
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5 Summary and Discussion 



Neutrinoless double beta decay is a crucial observable in search for physics beyond the 
Standard Model as it tests the fundamental symmetry of lepton number. The violation 
of lepton number is predicted in many models of new physics and most prominently, it 
provides the only probe of the absolute mass scale of light Majorana neutrinos. In this 
context, searches for Ou/Sp are highly complementary to neutrino oscillation experiments, 
direct neutrino mass determinations in Tritium decay and cosmological observations of the 
impact of neutrinos on large scale structure formation. If the recent results at the LHC 
pointing to a Higgs boson with a mass of about 125 GeV are confirmed, we can expect to 
make a giant leap in understanding the nature of mass generation at the electroweak scale. 
Even then, the nature and the lightness of neutrinos would still be unexplained and so far 
neutrinoless double beta decay is the only realistic probe to distinguish between the Dirac 
or Majorana nature of light neutrinos. 

As outlined in Section |3l there is large number of possible effective operators that can 
give rise to OuPP decay. Only a small selection of fundamental physics models that give 
rise to such LNV operators were presented in Section |H Given these ambiguities, a crucial 
problem is to distinguish between different mechanisms. One possibility is to compare 
results from Oz//3/3 with other neutrino experiments and cosmological observations. More 
generally, searches for physics beyond the Standard Model at other experiments, such as 
the LHC, can be correlated with Oz//3/3 in specific models. More directly, it is also possible 
to infer the dominant mechanism purely within the context of neutrinoless double beta 
decay and closely related processes. Relevant techniques discussed in the literature include 
(i) the comparison of Oz//3/3 decay rates in different isotopes |104H106j (this possibility is 
discussed in detail in the contribution by Fogli et al.), (ii) the comparison of Oz//3~/3~ 
with OZ//3+/3"'" |107] . (iii) the comparison of Oz//3/3 with electron capture |107] . (iv) the 
comparison of of Oz//3/3 decay to the ground state and an excited state [108] and (v) the 
experimental determination of the angular and energy distribution of the outgoing electrons 

in OZ//3/3 nmnMH]. 

Here we highlight the potential of measuring the angular and energy distribution 
in OuPP tracking experiments such as SuperNEMO |113j . Both the angular and energy 
correlation of the electrons depends on the effective operator mediating neutrinoless double 
beta decay. In the standard mass mechanism with V — A couplings, the electrons are 
mostly emitted back to back with comparable energies. This is very different from the 
^v+A process shown in Fig. H] (left), where the electrons are preferably emitted in the same 
direction with one electron taking away the majority of the energy. Fig. shows the 
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Figure 10: Simulated constraints at one standard deviation on the effective neutrino mass 
nil, and LRSM coupling (A) = ^v+a ^^r the isotope ^^Se at SuperNEMO from: (1) an 
observation of Oz//3/3 decay half-life at T1/2 = 10^^ y (outer blue elliptical contour) and 10^^ y 
(inner blue elliptical contour); (2) reconstruction of the angular (outer, lighter green) and 
energy difference (inner, darker green) distribution shape; (3) combination of (1) and (2) 
(red contours). The strength of the A contribution is assumed to be 30% of the standard 
mass mechanism. NME uncertainties are assumed to be 30% and experimental statistical 
uncertainties are determined from the simulation in [ 113] . 

impact of combining the total Ou/SP rate measurement with a determination of the angular 
and energy distribution shape at SuperNEMO. In the hypothetical scenario considered, 
both the standard mass mechanism and the operator of Cy^^ are assumed to contribute; 
whereas the total rate alone can not differentiate between the two contributions, taking 
into account the angular and energy shape can allow to pinpoint their relative size. 

If neutrinoless double beta decay is observed in next generation experiments at a 
level corresponding to an effective OuPP mass rriee ^ 10~^ eV, it would be an indication 
that the decay is caused by the exchange of light neutrinos, especially if corroborated by 
cosmological observations. On the other hand, such a conclusion is not straightforward as 
there is a large number of models which can trigger neutrinoless double beta decay. As 
discussed in this article, such new physics mechanisms can be economically categorized in 
terms of effective Lorentz invariant operators. Due to the black box theorem [2l[3l[T3], the 
observation of neutrinoless double beta decay will prove that the light left-handed neutrinos 
are Majorana particles, but from this alone it is not possible to infer the dominant operator 
and mechanism. 
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